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Abstract 

The theory of perfect fluids is reconsidered from the point of view of a 
covariant Lagrangian theory. It has been shown that the Euler-Lagrange 
equations for a perfect fluid could be found in spaces with afflne connec- 
tions and metrics from an unconstrained variational principle by the use 
of the method of Lagrangians with covariant derivatives (MLCD) and ad- 
ditional conditions for reparametrizations of the proper time of the mass 
elements (particles) of the perfect fluid. The last conditions are not related 
to the variational principle and are not considered as constraints used in 
the process of variations. The application of the whole structure of a La- 
grangian theory with an appropriate choice of a Lagrangian invariant as 
the pressure of the fluid shows that the Euler-Lagrange equations with their 
corresponding energy-momentum tensors lead to Navier- Stokes' equation 
identical with the Euler equation for a perfect fluid in a space with one 
afflne connection and metrics. The Navier-Stokes equations appear as 
higher order equations with respect to the Euler-Lagrange equations. 



1 Introduction 



1. Every classical field theory could be considered as a theory of a continuous 
media £Q 4- 0] with its kinematic characteristics. On the one side, a classical 
field theory is usually based on three essential structures: 

(a) The Lagrangian density, 

(b) The Euler-Lagrangian equations, 

(c) The energy-momentum tensors. 

The structures in (a) - (c) determine the scheme of a Lagrangian formalism 
for describing the properties and the evolution of a dynamic system. The find- 
ing out of the Euler-Lagrange equations and the energy-momentum tensors is 
related to a variational principle applied under certain mathematical conditions 
for the variables and their variations [S], @|. 
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On the other side, a theory of a fluid is based on 

(a) A state equation 

(b) A variational principle with a given Lagrangian density depending on 
the variables constructing the equation of state (rest mass density, entropy, 
enthalpy, temperature, pressure etc.) and the velocity of the points of the fluid. 

2. Since a theory of a fluid is a special case of a theory of continuous media, it 
could be worked out by means of the same mathematical tools used in a classical 
field theory. Moreover, the Lagrangian invariant L in the construction of the 
Lagrangian density L := \J—d g ■ L, d g := det(gij) < 0, could be interpreted 
as the pressure p of a physical system, described by its Lagrangian density L. 
Here, the tensor g = g^ ■ dx l .dx j , dx l .dx^ := | • {dx l ® dx j + dx j (g) dx l ) is the 
metric tensor field in a diffcrentiable manifold M, dimM — n, considered as a 
model of the space {n — 3) or of the space-time (n = 4). The fact that the 
pressure p could be used as Lagrangian invariant L — p has been observed by 
many authors [B] (see the citations there) but used only for very special forms 
of the pressure p. The application of the assumption L :— p could be extended 
to many classical field theories and to theories of different types of fluids @] . 

3. Every Lagrangian density depends on a set of field variables. Part of 
them (or all of them) are considered as dynamic field variables, i.e. as variables 
which variations lead to the Euler-Lagrange equations or to energy-momentum 
tensors. If a field variable is not varied it is called a non-dynamic field vari- 
able. It depends on the problem to be solved which variables are considered as 
dynamic and which as non-dynamic field variables. Only in the case of find- 
ing out the energy-momentum tensors the Lie variation j^j of all field variables 
should be taken into account. If all dynamic field variables are independent to 
each other, there are no additional conditions between them. The variables are 
varied without any constraints and the variational principle used in that case is 
called unconstrained variational principle. If a part of the dynamic field vari- 
ables are depending in some way to each other there are relations (constraints) 
between them. The variation of the variables have to include the constraints 
between them and the variational principle used in this case is called constrained 
variational principle. In many cases the constraints have no obvious physical 
interpretation or lead to additional problems. There are attempts for finding 
out conditions which could show us which field theory could be derived from an 
unconstrained variational principles and which cannot [§] ■ On this basis Schutz 
and Sorkin came to the conclusion that the theory of a perfect fluid could not 
be derived from an unconstrained variational principle. At that, they used the 
method of Lagrangians with partial derivatives (MLPD) [5]. The theory of a 
perfect fluid [related to the Euler equation in hydrodynamics (1727-1741)] is 
considered as the oldest of classical field theories |S] and it is worth to be inves- 
tigated by the use of some of the modern mathematical tools with respect to its 
structures and elements. 

4. In the present paper we have used the method of Lagrangians with co- 
variant derivatives (MLCD) instead of the method of Lagrangians with partial 
derivatives (MLPD) to reconsider the problem of derivation of the theory of a 
perfect fluid from an unconstrained variational principle. It has been shown that 



2 



the Euler-Lagrange equations for a perfect fluid could be found in spaces with 
affine connections and metrics 3 from an unconstrained variational principle by 
the use of the method of Lagrangians with covariant derivatives (MLCD) and 
additional conditions for reparametrizations of the proper time of the mass ele- 
ments (particles) of the perfect fluid. The last conditions are not related to the 
variational principle and are not considered as constraints used in the process of 
variations. The application of the whole structure of a Lagrangian theory with 
an appropriate choice of a Lagrangian invariant as the pressure of the fluid shows 
that the Euler-Lagrange equations with their corresponding energy-momentum 
tensors lead to Navier-Stokes' equations containing the Euler equation for a 
perfect fluid as a special case in a space with afnne connections and metrics [7], 
3 . The Navier-Stokes equation appear as higher order equations with respect 
to the Euler-Lagrange equations. 

5. The structure of the paper is as follows. In Section 2 a brief recall of 
the method of Lagrangians with covariant derivatives (MLCD) is made and ap- 
plied to a Lagrangian density depending on a metric, a scalar field, and on a 
contravariant non-isotropic (non-null) vector field. In Section 3 the invariant 
projections of the energy-momentum tensors are considered for the correspond- 
ing Lagrangian density. In Section 4 the covariant divergencies of the energy- 
momentum tensors and their relations to Navier-Stokes' equations are given. In 
Section 5 the invariance of dynamic characteristics under changing the proper 
time of the trajectory of a mass element (particle) in a fluid is considered. In 
Section 6 the notion of perfect fluid and its mathematical models arc considered. 
In Section 7 the method of Lagrangians with covariant derivatives is applied to 
a special type of a Lagrangian density related to a perfect fluid in special types 
of (L n , <?)-spaces and in (L n , g)-spaces and the corresponding Euler-Lagrange 
equations and energy-momentum tensors are found. In Section 8 the Navier- 
Stokes equation is considered in (L n , (^-spaces. It is shown that this equation 
is identical with the Euler equation if considered for a perfect fluid with the 
given Lagrangian density in a (L n , g)-space. The final Section 9 comprises some 
concluding remarks. The most considerations are given in details (even in full 
details) for these readers who are not familiar with the considered problems. 
If some details or abbreviations are not defined or introduced in the paper the 
reader is kindly asked to refer to |2j and 0] where all preliminary informations 
for reading this paper are given. 

2 Method of Lagrangians with covariant deriva- 
tives (MLCD) 

2.1 General remarks 

The method of Lagrangians with covariant derivatives (MLCD) is considered in 
details in 0, jH], 0] We will recall here only the main structures and properties 
of this method. 
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The method of Lagrangians with covariant derivatives (MLC'D) is a La- 
grangian formalism for tensor fields based on: 

(a) A Lagrangian density L of type 

L = -d g .L(gij, g i:j . k , g ij;k[ i,V A B ,V A B^ yA B-.i-.j) > W 

where g^ are the components of the covariant metric tensor field g, denotes 
covariant derivative with respect to the co-ordinate x k (or with respect to the 
basic co-ordinate vector field dk), V A b are components of tensor fields V € 
<X> fc i(M) with finite rank, A and B are co-operative indices. 

(b) The action 5 of a Lagrangian system described by means of the La- 
grangian density L 

S = j L.d^x = J L.duj , (2) 

v„ v„ 

where V n is a volume in the manifold M with dimM = n; dui = ^—d g .d^ n 'x is 
the invariant volume element. 

(c) The functional variation 6S of the action S with the condition for the 
existence of an extremum 

SS = S j L.d^x = j SL.d^x = . (3) 

v„ v n 

(d) The functional variation of the Lagrangian density L in the form 

dL = - — .6 gij + ■o{gij;k) + „ ■<HS»j;/w)+ 

ar fir ffr 

(e) The variational operator 5 obeying the commutation relations leading to 
commutation of S with the covariant derivatives: 



S(gij-k) = (5gij);k , 8(gij;k-,i) = {8gij);k;i , 
5{V A B-,i) = (6V A B ) ;i , S(V A B .i.j) = (6V A B ). Aj . 



(5) 



(f) The Lie variation £^S of the action S 

£^S = £t J L.d^x = J I^L.d^x . (6) 

(g) The Lie variation of the Lagrangian density L in the form 



4 



Remark. In the MLCD [because of the commutation relations] the affine 
connections T and P, and their corresponding curvature tensors are considered 
as non-dynamic field variables (ST^ k = 0, SP^ k = 0, SR Z jki = 0, 5P l ju — 0). 
Therefore, in the MLCD variations of the components and their covariant deriva- 
tives of the covariant metric tensor g and the non-metric tensor fields V are 
considered for given (fixed) affine connections and for fixed and determined by 
them types of transports of the tensor fields. Of course, the question arises how 
the affine connections can be found if not by means of a Lagrangian formalism. 
The first simple answer is: the affine connections (or the equations for them as 
functions of the co-ordinates in M) can be found on the grounds of the method 
of Lagrangians with partial derivatives (MLPD) and then the components of 
the tensor fields (as functions of the co-ordinates in M) can be determined by 
means of MLCD. This answer could induce an other question: why two meth- 
ods have to be applied when one is enough for finding out all equations for all 
dynamic field variables. There are at least two possible answers to this question: 
1. The MLCD ensure the finding out equations for tensor fields (as dynamic 
field variables). These equations are (a) covariant (tensorial) equations and (b) 
form-invariant (gauge invariant) equations with respect to the affine connec- 
tions. The affine connections, related to the type of the frame of reference [5] 
could be determined on the grounds of additional conditions and not exactly by 
means of a variational principle. 2. The MLPD can ensure the consideration of 
the affine connections as dynamic field variables and the finding out their field 
equations. It cannot give direct answer for the type (tensorial or non-tensorial) 
of the equations obtained for the tensor field variables and their corresponding 
energy- momentum quantities. The tensorial character of quantities and rela- 
tions concerning the tensor fields variables has to be proved (which, in general, 
could be a matter of some difficulty) |1U| . 

2.2 Method of Lagrangians with covariant derivatives ap- 
plied to a scalar and to a vector field 

2.2.1 Lagrangian density 

1. The mathematical description of a state of a fluid is usually made by means 
of the velocity vector field u of the moving mass elements (particles) of the 
fluid and other two functions (for instance, the pressure p and the rest mass 
density p) If u, p, and p are given then the state of a fluid is determined. 
This fact gives rise to the idea a Lagrangian invariant L = p depending on the 
velocity vector field u and on the rest mass density p to be used in the MLCD 
for finding out the Euler-Lagrange equations and their corresponding energy- 
momentum tensors in a space with affine connections and metrics [(L n , g)-space] 
(considered as a model of space or space-time). 
2. Let a Lagrangian density of the type 

L := yj~d~ g ■ p{g l} , p, p ;i = p A , u l ) (8) 

be given, where d g := det(gij) < 0, gij = gji(x k ) are the components of a covari- 
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ant metric g, p is an invariant function on M, p £ C r (M), r > 2, p £ <g)° o(M), 
p. i = p i := dp/dx l , u l are the components of a non-null (non-isotropic) vector 
field u £ T(M) in a co-ordinate basis {<%} C T(M), u = u l -di, g(u, u) := e 0, 
i,j = 1, n, dimM = n. The Lagrangian invariant L := p could be interpreted 
as the pressure p of the fluid. The invariant function p could be interpreted as 
the rest mass density of the fluid corresponding to one of its energy-momentum 
tensors 5 . The contravariant non-isotropic (non-null) vector held u and its 
components u l could be interpreted as the velocity vector field of the mass 
elements (particles) of the fluid. 

By the use of the MLCD we can find the Euler-Lagrange equations and the 
energy-momentum tensors corresponding to L. The Lagrangian density L is a 
degenerated Lagrangian density with respect to u l and gtj, i.e. L contains no 
derivatives of u and gij with respect to the co-ordinates x % or to the basic vector 
fields a = di. 



2.2.2 Euler-Lagrange's equations 

(a) Euler-Lagrange's equations for the rest mass density p 

The ELE for the rest mass density p could be found in the form jS], 0] 



Sp 



p = Q , (9) 



where 



P = di-p- , qi = T ki k -l-g kl -g k i;i+g l k -gf ;i ■ (11) 

Therefore, the ELE for p could be written as 
dp , dp dp 

Tp-^+^-jp-r ■ (12) 

(b ) Euler-Lagrange 's equations for the vector field u 

The Lagrangian density L is degenerated with respect to the vector field u, 
i.e. the ELEs could be written in the form 

(c) Euler-Lagrange's equations for the metric tensor g 

The Lagrangian density L is degenerated with respect to the components 
gij of the tensor field g. Therefore, 

S g P _ dp 
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and wc have the ELEs in the form 



" P +Lp.gV=0 . (14) 



dgij 2 

If the ELEs for gij are fulfilled then the pressure p should obey the equation 

dp n 2 dp 

9ij + 7i-P = : P = o 9ij • (15) 



dg^ 2 ' n dg i3 

The last condition means that p should be a homogeneous function of gij 
of degree m = — (2/n). If the pressure p does not have this property then the 
metric tensor g could not be considered as a dynamic variable, i.e. it should not 
be varied. 

2.2.3 Energy-momentum tensors 

In accordance to the general scheme of the MLCD we could find the explicit 
form of the corresponding to the pressure p energy-momentum tensors. 

(a) Generalized canonical energy-momentum tensor 

<>. J i^ -i>:, /-//: • (16) 

(b) Symmetric energy-momentum tensor of Belinfante 

sTi^-p-gi . (17) 

(c) Variational energy-momentum tensor of Euler- Lagrange 



7n j &uP k j I $gP J m S gP i m 

dp j dp j m dp j m 

= ft?-* 1 -^^^ -9ml- ^-91-9-, - 9km ■ 

Because of g Jk = f m i ■ g mk = f m j • g km = g kl , g kl = g lk , gf ■ g ml = 
/" i ■ 9n ■ 9ml —~ f i 9nl = 9ji we have 

7T i d P j o d P j noi 

3 Invariant projections of the energy- momentum 
tensors 

Since the vector field u could be interpreted as the velocity of the mass elements 
(particles) of a fluid, we can project the energy- momentum tensors 6i 3 , s Ti 3 , 
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and Q i J by the use of the vector field u and its corresponding covariant and 
contravariant projective metrics h u and h u 

1 1 

K = g -g(u)<gig(u) , h u = g u®u , (19) 

e e 



j > 



e = g(u,u):^0 , g = g v ■ d l .d. 
di.dj = - ■ (di <g> dj + dj <g> d l ) , 

g( u ) = 9i] ■ uJ ■ dx l = g {j ■ u J ■ dx l . 

By the use of u, h u , and h u we can find the different invariant projections 
of Qi^ , sTi? , and which have physical interpretation. 

3.1 Invariant projections of the generalized canonical energy- 
momentum tensor 9 

1. Rest mass density p g . 

Pe = ^ • Ok 1 ■ uj ■ u k = 

1 k / dp iN 

= e2' Ul ' U '^dp~.' P;k-P-9k) = 

= ^•^•"-•^•P; fe -P-^- 1 )- (20) 
2. Conductive momentum density e W. 

7T -=~-k&i -uj-h = - ■ — • p.y h . (21) 



3. Conductive energy flux density e • s 

1 uH „ .a.k J _ 1 U ij „ _&P 

dp 



W ■ g Tk - kBi k -u l = - ■ g Tk - ^ ■ p,-u l . (22) 



>■"« h ' J -"--^-^-^ ■ (23) 
4. Stress tensor density e S. 

^ i :=h*- 9rrk 0j-h™* = h* k -g m -^-p. m -h™* . (24) 

3.2 Invariant projections of the symmetric energy-momentum 
tensor of Belinfante 

1. Rest mass density p T . 

Pt~ ^' sT k i -uj-u k = --L..p. M _. u fc . (25) 
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2. Conductive momentum density n. 

T W i := - ■ Ti k -u-j:-h u = , Ti k = . (26) 

3. Conductive energy flux density e • T s. 

T s l := - • h y ' • • Ti k ■ u l = , T ; fc - . (27) 

e J 

4. Stress tensor density T S. 

T S ij :=h ik - m - T m l -h mi = . (28) 

3.3 Invariant projections of the variational energy-momentum 
tensor of Euler-Lagrange 

1. Rest mass density pq. 

1 7^ i k 2 d P i k 

PQ = —Z-Qk - U r U = "2 ' «— ' ffkrn ' 9n ' u i ' 11 ~ 

e 2 <9w fc 1 

2 dp t 1 dp t 

^ = ?-%^-^- U "- U --e^- U • (29) 

2. Conductive momentum density ®ir. 

3. Conductive energy flux density e ■ ®s. 
Q s l = ---hv-ffj,;- Qt k -u l = 



e ^ 



1 i a J r dp fen dp k 

2 . ™ T k dp 1 , • J k d P 

Qs< ='-^'-^-^-^- 7 • (3D 
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4. Stress tensor density 

Qgij = _ h ik . Qj . 



= -^■m^ l -^-h mj + 2-h ik - m --^--g mr - g ^-h^ 

«S« =2-h ik - m --^--grnn- . 

dg n i 

3.4 Noether's identities 

From the second covariant Noether identity 

9- S T = Q 

we obtain the identities: 
!• Pe = Pt - Pq 

dp I dp t n dp t 

^■«T« -P^e.— k .u -2.— n9 - km .u n .u 

2. e W= T W-Qir: 

dp ,u _ dp - u dp - H 



dp, k 11 k du> dg, 

3. e s= T s-Qs: 



dp — dp — 



4. e S= T S~ ®S: 

hik -m-^-p-, m -h mj ^2.h ik . m .-^-. 9mn .h^ . 

3.5 Explicit form of the energy-momentum tensors 

1. Generalized canonical energy-momentum tensor 9. 

9 = (p g + i • k ■ L) ■ u <g> g(u) -p- Kr + u® g( e n) + 8 s <g> + ( e #)g . 

2. Symmetric energy-momentum tensor of Bclinfante S T. 

S T = — p • Kr . 
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3. Variational energy- momentum tensor of Euler-Lagrange Q. 

Q = -p Q -u®g{u)-u®g{ Q -rt)- Q s ®g{u)-( Q S)g . (40) 

The form of the symmetric energy-momentum tensor of Belinfante shows 
that a fluid, described by a Lagrangian density of the type L = \J—d g ■ p{gij, 
p, p.^ u l ), could be considered as a fluid in equilibrium state |12) . Such an 
equilibrium state could be related to the Euler equation of a static state of a 
fluid. 

4 Covariant divergencies of the energy-momentum 
tensors and Navier-Stokes' equations 

4.1 Covariant divergency 

The covariant divergency of the energy-momentum tensors could be written in 
the forms respectively 0, 0] 

1. Covariant divergency of the generalized canonical energy- momentum ten- 
sor. 

g(S§) = (p e + - -k-p)-a + 

e 

+ \u(p g + - ■ k ■ p) + (p g + - ■ k ■ p) ■ 5u + S e s] ■ u — 

e e 

-g(Krp) - p ■ g(6Kr) + Su ■ e W+ V u e ir + We^u + 

+ {p e + - ■ k ■ p) -g(W u g)(u) +g(W u g)( e W) +g(Wo^g)(u) + 

e 

+9(S(( e S)g)) . (41) 

2. Covariant divergency of the symmetric energy-momentum tensor of Be- 
linfante. 

5 S T = —5(p-Kr) — ~(Krp + p-SKr) , 
g(S s T) - -g{Krp)-p-g{5Kr) . (42) 

3. Covariant divergency of the variational energy-momentum tensor of Euler- 
Lagrange. 

g(SQ) = -p Q ■ a- (up Q + p Q ■ Su + S Q s) ■ u- 
-Su ■ q tt - V„ q tt - Vq s « - 

-Pq ■ 9{Vu9){u) - g{V u g)^T:) - g(W Qs g)(u) - g(5(( Q S)g)) (43) 

4.2 Navier-Stokes' equations 

The Navier-Stokes equations follow from the second Navier-Stokes identity @] 

hu[g(F)] + h u [g{Se)]=0 (44) 
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in the form 

h u [g(S9)}=0 . (45) 
In its explicit form the Navier-Stokes equation could be written as 

{pe + ^-k-p}- h u (a) - h u [g(Krp)} - p ■ h u [g(5Kr)] + 

+5u ■ h u ( e W) + h u (V u e Jf) + h u {V<,- s u) + 

+ {pe + \-k-p)- h u [g{V u g){u)] + h u [g(\7 u g)( e w)} + 

+h u [g(Ve ig )(u)} + h u [g(S(( e S)g))} =0 . (46) 
If we write the explicit form of the expressions 

K(a) = g(a) - i • g(u,a) ■ g(u) , (47) 

h u [g(Krp)] = g[g(Krp)} - i • g{u,g{Krp)) ■ g{u) = 

= Krp - - • g(u, g(Krp)) ■ g{u) (48) 
e 

in the Navier-Stokes equation we obtain it in the form 

(Pe + ~ ■ k ■ P) ■ 9(a) ~ ~ ■ (Pe + ~ ■ k ■ P) ■ 9(u, a) ■ g{u) - 
e e e 

-Krp+ - ■ g(u,g(Krp)) ■ g(u) -p- h u [g(5Kr)} + 

e 

+5u ■ h u ( e n) + h u {V u 9 7f) + huiVe^u) + 

+{p e + - ■ k -p) ■ h u [g{V u g)(u)] + h u [g(W u g)( e w)} + 



+h u [g(Ve^g)(u)} + h u [g(6(( e S)g))} =0 . (49) 

After contracting the last (previous) equation with g and introducing the 
abbreviation 

p:=(p e + --k-p) (50) 

e 

we find the Navier-Stokes equation in the form 

p ■ a - g(Krp) - - ■ [p ■ g(u, a) - g(u,g(Krp))] ■ u - 

-P ■ g(h u )[g(5Kr)} + Su ■ g(h u )( W) + g(h u )(V u e W) + g{h u ){Ve- s u) + 

+p ■ g{K)[g{V u g){u)] + g(h u )[g(V u g)( e if)}+ 



+g(h u )[g(Vo^g)(u)}+g(h u )[g(S(( e S)g))} =0 . (51) 
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Special case: (L n , c/)-spaces. e n := 0, e s := 0, e S := 0. 

p ■ a - g(Krp) - ^ • [p ■ g(u, a) - g(u, g(Krp))} ■ u- 

-p-g(h u )[g(SKr)}+p-g(h u )[g(V u g)(u)} = . (52) 
Since g{h u )g = h u , we have 

p-a- g(Krp) - ^ • [p ■ g(u, a) - g(u, g(Krp))} ■ u- 

-p ■ h u (SKr) + p ■ h u [(V u g)(u)} = . (53) 

Special case: U n - and y„-spaces. V^g = for e T(M), 6 Tf :— 0, e s := 0, 
"S: II. 

p-a- g(Krp) --■[?■ 9(u, a) - g{u, g(Krp))} ■ u- 

-p-h u {5Kr) =0 . (54) 

Special case: U n - and F„-spaces. V^g = for V£ G T(M), 6 tt := 0, e s := 0, 
e S := 0. 

Since V M g = and SKr = 0we have the Navier-Stokes equation in the form 

p-a- g{Krp) - - ■ [p ■ g(u, a) - g(u, g(Krp))} -u = . (55) 
e 

The contravariant non-null (non-isotropic) vector field u is interpreted as the 
velocity vector field of the moving mass elements (material points, particles) in 
the fluid. The vector field a is interpreted as the acceleration of the same mass 
elements of the fluid. If we write the Navier-Stokes equation in a basis (in a 
co-ordinate or non-co-ordinate basis) we obtain its form 

p-a 1 - g t] ■ p d - - • [p • g- r u k ■ a 1 - u l ■ pj] ■ u l = . (56) 
On the other side, 

g(u,a) = - ■ [ue - (V u g){u,u)] , a = X7 u u , 
V„e = V u \g(u,u)]=ue = (V u g){u,u) + 2-g(a,u) . 

4.2.1 Navier-Stokes' equation in U n - and V„-spaces 

In a U n - or V^-space we obtain the Navier-Stokes equation in the form 

11 

p ■ a 1 - g v ■ pj [- • p ■ ej ■ u° - v? ■ p tj ] ■ u 1 = . (57) 
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If the last term of the type of / • u vanishes then the expression takes the 
form of the Euler equation 

T>-"' •/'..- ■ (58) 

The vanishing of the term / • u with 

/ := - • [- • p ■ ej -u 3 -u 3 ■ p d ] (59) 

could be achieved in two different ways: 

(a) By additional conditions (constraints) on the vector field u and on the 
pressure p. If e := const. ^ and J 2 := (for u = 4-) then / = 0. These 
additional conditions should be included as constraints on the Lagrangian in- 
variant p or introduced by definition of the pressure p and the velocity u. The 
condition e := const. ^ is usually assumed in the relativistic mechanics in 
V^-spaces. = should be in some way included in the Lagrangian formalism 
for p. This is the reason for the statement that the Euler equation could not be 
found on the basis of an unconstrained variational principle. 

(b) By changing the proper time r as the parameter of the trajectory (world 
line) x l (t) of a mass element. This change is not related to a variational principle 
and could be performed after all steps related to the Lagrangian formalism. The 
change of the proper time r is related to the canonical form of an auto-parallel 
(or geodesic in a V^-space) equation. 

Let us now consider the consequences from the changing the proper time 
r in an equation of motion of a mass element (particle) in the Newtonian and 
relativistic dynamics. 

5 Invariance of dynamic characteristics under 
changing the proper time of the trajectory of 
a mass element (particle) 

1. One interesting question is related to the invariance of physical laws under 
time's transformations. The general covariance of a dynamic law under co- 
ordinate transformations (diffeomorphisms of the differentiable manifolds con- 
sidered as models of space or space-time) is well known and used in different 
physical theories. Let us now consider how a dynamic law changes under a 
change of the time's parameter in it. A good example for the behavior of an 
equation under changing the time's parameter in it is the auto-parallel equa- 
tion in spaces with affine connections and metrics. It is well known that the 
change of the proper time is related to the canonical or non-canonical form of 
an auto-parallel equation ^3], [3], 0]. Nevertheless, the physical interpreta- 
tion of a change of the proper time is not considered in details. What is the 
physical interpretation of the change of the time's parameter in Newtonian and 
relativistic physics? 
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5.1 Transformation of the time's parameter in the second 
law of the Newtonian mechanics 

The second law in the Newton mechanics could be written in the form 

V u p = F , . P = p-u , p = p(x k ,u,t) , F = F(x k , u, t) . 

(60) 

If the parameter t interpreted as time's parameter is changed with a new 
parameter r with 

t = t(r), T = r(t), t,TeC r (M), r>2 , dimM = n (n = 3J6}) 
^ = a(x k ) = a(x k (r))=a(x k (t)) ,u:=±=a-u=^-j t , (62) 

F = F(x k ,u,t) = F(x k ,- -u , t( r )) . (63) 

a 

Then the second law, written by the use of the new time's parameter r will 
have the form 

V u (p- u) = • \7u(p-u) - \ -p- (ua) ■ u = F(x k , - -u , t(r)) . (64) 



Therefore, 



V„p = • V^p - \ ■ p ■ (ua) ■ u = F , (65) 



or ar 

Vup = p ■ (ua) -u + a 2 ■ F = F + J -u , (66) 
VuP = J-u + F . (67) 

where 

p = p -u , / = p ■ (ua) , F = a 2 ■ F . 

The change of the time's parameter t in the Newton second law generates 
in general an additional force proportional to the velocity u (or u). A force of 
this type (proportional to the velocity) is usually interpreted as a friction force. 
Only under a transformation t = t(r) for which 

1 1 d . dt. 1 d 2 t 

ua = — ■ ua = — ■—( — ) = — ■ —r — , a^O , (68) 
a a dr dr a dr z 

the second Newton law remains invariant up to the constant function (1/a 2 ) 

V uP =-^-VvP = F . (69) 
cr 

The condition 

leading to the conditions 

= do = const., t — a • r + b a , b a = const. , (71) 
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shows that only under linear transformations of the time t = t(r) the Newton 
second law is form invariant 

V uP = F , V v p = F . (72) 

If a mass element (particle) with rest mass density p is a free test particle 
(F = 0) moving under condition 

V uP = , (73) 

then the change of the parameter t generates a friction force a ■ f ■ u and the 
considered motion of the particle is no more a free motion but a motion un- 
der a friction force proportional to its velocity. Therefore, we can generate or 
remove friction forces in Newtonian mechanics by performing an appropriate 
transformation of the time's parameter t. 

Remark. Since u = a ■ u, the Lie derivative of u along u is £ u u — [u, u] — 
(not) ■ u. The dragging of u along u should not change u if we consider a free 
motion of a particle. This is so, because of the equation of motion \7 u u = 
leading to V^ru = and to the condition V„u = 0. Therefore, the Lie derivative 
£ u u = V u u — Vu« in a E n - or V^-space (n — 3,4) should vanish. The Lie 
derivative £ u u = —£uU vanishes if and only if ua — or ua = 0. This 
means that the function a should be a constant function, different from zero, 
on trajectories with tangent vectors u and u. Under this condition (ua = 0) the 
momentum of a free moving particle is an invariant quantity V u p = VuP under 
changing the time's parameter of the trajectory of the particle. 

5.2 Transformation of the time's parameter for a free mov- 
ing spinless test particle in spaces with affine connec- 
tions and metrics 

1. If we consider the auto-parallel equation as the equation for description of the 
motion of a free spinless mass element (test particle) with constant rest mass 
density we should distinguish its two forms [13], ^1] 

(a) Canonical form 

V u u = , (74) 

(b) Non-canonical form 

V u u = f-u . (75) 

2. The transition from the canonical to the non-canonical form could be 
realized to a transformation of the proper time (parameter) of the trajectory of 
the particle. 

If u := d/dr and u := (I /a) ■ u := d/df then 

V u « = \7 a .u(a -u) — a ■ [(ua) ■ u + a ■ V^w] = 

= a 2 ■ Vu& + (ua) ■ u = a 2 - a + / • u , (76) 
a = f ■ u + a 2 - a , / = ua , a — Vutt . (77) 
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The acceleration a after changing the parameter t is expressed by two terms: 
/ • u (proportional to u) and a 2 ■ a proportional to the acceleration a defined by 
the new proper time f. The term / • u defines a friction acceleration. 

3. If the auto-parallel equation is given in its non-canonical form as V u u = 
f ■ u with u = d/dr then we can change the parameter r with a new parameter 
f := t(t), r = t(t), so that 

_ d 1 dr d dr 1 

df a df dr df a 

The auto-parallel equation V u u = / • u could be written in the forms 

V„u = Va-n(a -u) = f -a-u , (79) 

[a ■ (ua) - / • a] ■ u + a 2 ■ V^u = . (80) 
The function a could be chosen in such a way that the condition 

a ■ (ua) - J ■ a = (81) 

could be fulfilled. Since a :^ 0, we have 

ua — f = Q:a = ao+ff-df , «o = const., (82) 
7 = 7(r)=7(r(r)) = /(r) , 



or 

a = So • exp[ J f(r) ■ dr] , So = const. . (83) 

If a fulfills the previous condition then V„t( = and the auto-parallel equa- 
tion is represented in its canonical form. 

Therefore, an acceleration of the type of a friction acceleration (proportional 
to the velocity vector field of a mass element) could be generated or removed by 
the use of transformation of the proper time of the mass element. This fact could 
be used for removing terms proportional to the velocity in the structures of the 
energy-momentum tensors and in the structure of Navier-Stokes' equation. 

4. Let as now assume that the acceleration a of a mass element with rest 
mass density p fulfills the equation 

V u u = a = f-u , f £C r (M) , r > 2 . 

Then the change of the momentum p = p ■ u of a mass element moving with 
the above acceleration could be written in the form 

V„p = V u (p ■ u) = (up) ■ u + p ■ a = (84) 
= (up) ■ u + p ■ f ■ u = (up + f ■ p) ■ u . 
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If the momentum p does not change under the acceleration a = f ■ u, i.e. if 

V M p- V u (p-u) = (85) 

then 

up + f-p = , u^O , p^O , (86) 
and the rest mass density p should obey the equation 

~~ = ~f ' U ( 1 °9P) = ~f . d ( l °9p) = -f-dr , 

and therefore, 

P = Po ■ ex P[- J f( T ) ■ dT \ , Po = const - ( 87 ) 
Special case: If / := A := const. :> 0, then 

P = Po ■ ^ X ' T ■ (88) 

5.2.1 Law of self-decay of a system of n-particles in a (L n , g)-space 

Let us now assume as above that the rest mass density p of a mass element with 
equation of motion of the type V u u = / • u is a composition of the constant rest 
mass density p of n particles, i.e. 



p := n ■ p , p = const. (89) 



Then 



V u (p-u) = V„(n- p -u)= p -(un)-u+ p -n-V u u = 
= ( un ) ■ Po ' u + Po ■ n ■ f ■ u = 

= [un + n ■ /] • po • u . (90) 

If V u p = V„(p ■ n ■ u) = 0, i.e. if the momentum of the mass element with 
momentum p and with rest mass density p = n ■ p docs not change along its 
trajectory (world line), then 

un + n ■ f = (91) 

and 

n = no ■ exp\— I f ■ dr] , n = const. (92) 

The number n is the number of the particles in the mass element p at the 
time t. no is the number of the particles in the mass element p at the time r . 

n{r)=n{T )-e l - J -o f - d ^ . (93) 
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Special case: f := A := const. > 0. 

n(r) = n(r ) • e~ A ' T . (94) 

The last expression is exactly the law of self-decay of a system of n particles in 
the time r. Therefore, if a mass element with rest mass density p (consisting of 
n particles with rest mass density p — const.) does not change its momentum 
density p = p-u along its world line with tangent vector u under an acceleration 
of the friction type a — f ■ u with / = A = const, then the number n of the 
particles in the mass element should obey the law of self-decay 

n(r) = n(r ) • e~ A ' T , n — const. , A = const. (95) 

It follows that from the one side, accelerations of the type a = f - u could be 
brought to zero by a corresponding transformation of the proper time r and on 
the other side, it could be important for the description of a self-decay law on 
a classical level. 

The consideration of the law of self-decay on an elementary mechanical level 
shows that the reparamctrization of the proper time could have important for 
the classical physics corollaries. The motion under friction acceleration and 
preservation of the momentum of a mass element despite of this friction accel- 
eration (which can be removed by a proper time reparametrization) could lead 
on the one side, to a self-decay of a mass element (emission of particles), but 
on the other side, it could lead to absorption of particles from the vicinity of a 
mass element (/ := A := const. < 0). The number of the emitted or absorbed 
particles n could also depend on the proper time r (related to the mass element 
p) in a more complicated manner than on the exponential law. 

5.3 Transformation of the proper time of a mass element 
in the Navier-Stokes equation in (L n , g)-spaces 

1. Now we can again consider the Navier-Stokes equation in the special case 
when % := 0, e s := 0, and e S := 0. It could be written in the form 

p ■ a - g(Krp) -J-u-p- h u (6Kr) + p ■ h u [(V u g)(u)} = , (96) 

or in the form 

p-a-g(Krp)-J-u-h u \p-6Kr-p-(V u g)(u)}=0 , (97) 

where 

/:= - e - [p-g(u,a)-g{u,g(Krp))} , (98) 

h u [p -SKr-p- (V u g)(u)} = p ■ h u (SKr) - p ■ h u [(V u g)(u)] . (99) 
The expression 

p-a-7-u = p-(a-l)=p-(a-f-u) , j = f - , (100) 
P P 
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could be transformed in an other expression by changing the time's parameter 
t, where u = d/dr and V„w = a, with a new time's parameter f — t(t). Then 

a — f ■ u = [a(ua) — f -a] -u + a 2 -a , a — V ¥ ti . (101) 

The function a (see above) could be chosen in such a way that the condition 

ua-f = (102) 

could be fulfilled. Then 

a = So • exp[ J /(r) • dr] , a = const., (103) 

and the Navier-Stokes equation could be written in the form 

p-a 2 -a-g(Krp)-h u [p-5Kr-p-(X7 u g)(u)] =0 , p :^ , a + . (104) 
or in the form 

a- = -^-g(Krp)- = -^-h u \p-6Kr-p-(V u g)(u)]=0 . (105) 
pa pa 

Special case: ?7„-and V^-spaces. S := C, V^.g := for V£ G T(M). For this 
type of spaces the Navier-Stokes equation has the form 

a = = ■ ^ • g(Krp) (106) 
p a 

because of SKr = and V u g = 0. 

2. The metric 7? is an arbitrary given metric. We can now introduce a 
conformal to g metric tensor g 

g~a 2 -g. (107) 
Then the corresponding contravariant metric g will have the form 

7j=a- 2 -g (108) 

because of the relations 

o? -g^-g a -g k3 = a 1 ■ g ji g- ! . 

g* = 7T ■■■■„>' . 

With the new introduced (given) conformal to g metric g the Navier-Stokes 
equation could be written as 

a=yTj{Krp) , (109) 
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or in a co-ordinate basis as 

=•//'•'•/'... ( 11Q ) 

which is exactly the Eulcr equation in (pseudo) Riemannian spaces with or 
without torsion. 

3. In a 3-dimensional Euclidean space E3 the Euler equation will have the 
form 

a b = yg bc -p, c , a,b,c= 1,2,3. (Ill) 

If we now chose the metric g as the flat metric in £"3, i.e. if g = r\ = 
t] ab ■ dx a .dx b , r] ab = (1, 1, 1) then we obtain the classical Euler equation 

a b = =-V bc -P,c , (112) 

a b = u b . c -u c , Tf:=^l = ^l , dr = dt, (113) 

af at 

which could be written in a form, usually used in the classical hydrodynamics 

u b A + u a -if ,a- = -p,b = $ , (H4) 
P 

_ h du b _ h du b 

u ' i = ^T ' u - = ' 

Therefore, the Navier-Stokes equation for a perfect fluid in a 3-dimensional 
Euclidean space of the classical hydrodynamics is identical with the Euler equa- 
tion for the perfect fluid. 

4. We can also write the Navier-Stokes equation in the form 

a-f-u-^-g(Krp)-h u [^-dKr-(W u g)(u)} = , (115) 

where ^ 

f=--[9{u,a)---g(u,g(Krp))] , 
e p 

p = p g + - -k- P . 

e 

If we change the parameter (proper time) r of the velocity vector field u = 
d/dr with a new parameter r = r(r), where r = t(t) with 

_ d dr d 1 

gJt df dr a 
dr 

u = a-u , a=— , aeC(M) , r > 1 , 
ar 



21 



then 

a = V u u = Va.u(a ■ u) — a 2 ■ a + (ua) ■ u , (H6) 

a = — i ■ [a — i • (ua) ■ u] , a^O . (H?) 

a a 

If a should be collinear to a (this means no change of the direction of the 
acceleration a if r is replaced with r) then 

a = ■ a and ua — . (118) 

a 2 

ua = means that the function a should not change along the trajectory 
with the parameter t. 

5. Let us now again consider the transformation of the expression 

a-j-u 

under the change of the proper time r. 

a — f ■ u = a 2 ■ a + (ua) ■ u — f ■ u = 
= a 2 ■ a + [(ua) — /]■« = 

= a 2 -a+ • (ua) - /] -u . (119) 
a 

The last term [i • (ua) — /] • u could vanish under two types of conditions: 

(a) Since a is until now an arbitrary given function a(r(x k )) = a(x k ) e 
C r (M), r > 1, we could specialize its form by the use of the condition 

ua- a- /= , a^O , (120) 

leading to the equation 

^ = / , u(loga) = .f , ±{log7S) = f , (121) 
with the solution 

a = a ■ exp( J f ■ dr) , a = const. (122) 

(b) If a should be collinear to a then ua = and a — / •u = a 2 -aif/ = 0. 
Remark. If a is collinear to u, i.e. a := fc • u then 

(fc — /) • u = a 2 ■ a + [i • (ua) — /] • u 
a 

and 

a 2 • a = [fc — / — — • (ua) + /] • u = [k — — ■ (ua)] ■ u 
a ' a 
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will be also collinear to u. Both equations appear as auto-parallel equations for 
the vector fields u and u collinear to each other. The auto-parallel equations 
could be considered as a very special type of Navier-Stokes' equation. They 
could be investigated separately from the more general cases when the trajecto- 
ries (world lines) of the mass elements are not auto-parallel trajectories for the 
velocities vectors u and u. 

Therefore, two possible conditions could be imposed for removing the term 
/ • u from the expression of the Navier-Stokes equation by the use of a transfor- 
mation of the proper time r 

(a) a — ao ■ exp(j f ■ dr) , «o = const. 

(b) ua = 0, / = 0. 

The condition ua = assures the collinearity of a and a even if / ^ 0. From 
physical point of view, condition (b) is more reasonable because of the condition 
ua = leading in any case to 



a, = ^jr ■ a or to a = a 2 ■ a . (123) 

a 



The condition / = has its physical consequences. It could be written in 
the form (for e :^ 0) 

g(u,a) = = -g(u,g(Krp)) . (124) 
P 

In a co-ordinate basis 

9iJ ■ u l ■ a 3 = i • gjj ■ u l ■ g jk ■ ■ p. m = 

= z-« ? -p,i ■ (125) 
P 

Special case: U n - and V^-spaces: S — C, \7^g = for V£ G T(M), e := 
g(u, u) := e = const. ^ 0. 

g(u,a) = ^■[ue-(V u g)(u,u)] = , (126) 
u'-pi - ^=0 . (127) 

If / = then the Navier-Stokes equation (when e n := 0, e s := 0, and 
°S := 0) takes the form of the generalized Euler equation in a (L„,g)-space 

a- = - g(Krp) -h u [^- SKr - (V u g){u)} - (128) 
P P 

or in a co-ordinate basis 

a^^-g^ -p^+hP ■{^■g k j , k -g jk ,i-u l -vF) . (129) 
P P 
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From the previous expression for a (or for a 1 ), it follows, after contraction 
with the vector field u, the condition 

g(u,a) = ^-g(u,g(Krp)) + [g(u)](h u )[^-6Kr-(V u g)(u)} . (130) 

Because of [g(u)](h u ) = 0we obtain 

g{u,a) = = -g(u,g{Krp)) . (131) 
P 

This means that the condition / = is automatically fulfilled if the gen- 
eralized Euler equation in a (L n , g)-space is valid. Therefore, the condition 
f = appears as a necessary and sufficient condition for the existence of the 
generalized Euler equation in a (L n , g)-space when 6 if :— 0, e s := 0, e S := 0. 

In conclusion, the Euler equation in a (L„, g)-space could be found as a 
special case of the Navier-Stokes equation under the following conditions 

e 7f = 0, 9 s = 0, e S = , (132) 
ua = , / = . (133) 

The last two conditions are related to the form invariancc of the acceleration a 
(within a scalar factor) with respect to the change of the proper time r of the 
world lines of the mass elements of a flow. At the same time these conditions 
appear as sufficient conditions for the existence of the Euler equation as a special 
case of the Navicr-Stokes equation if 7f := 0, e s := 0, e S := 0. 

5.4 Invariance of the projective metric h u corresponding 
to the contravariant non-isotropic (non-null) vector 
field u = d/dr under transformation of the time pa- 
rameter T 

1. If the time parameter r is transformed by the use of the transformation 

T = t(T),T =T(t), 

d dr d dr d 

u = — = — ■ — = a-u , a = — , u = — , 134 
dr dr df dr dr 

then the corresponding to u projective metrics 

K = 9 ~ ^ • g(u) ® g{u) , e = #(u,u):^0 , 
- 1 

n = q ■ u ® u , 

e 

do not change, i.e. h u and h u are form invariant under the transformation the 
time parameter r of the vector field u. On this basis we can prove the following 
propositions: 
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Proposition 1. The corresponding to a non- isotropic (non-null) contravariant 
vector field u = d/dr projective metrics h u and h u are invariant under the 
transformation of the parameter r by means of the transformation r = r(r) 
andr = r(r). 

Proof. From the explicit form of h u 

K = g-j ■ g(u)®g(u) 
and the form of u after the transformation of r 

u = a-u , t = t(t) , a :^ , 



it follows that 



/i„ = 5 ] • g(a ■ u) <g> g(a ■ u) 

g(a ■ u, a ■ u) 



1 

a 2 g(u,u) 



= 9 i^—:-g{u)®g{u) = hu . (135) 

g{u,u) 

Therefore, h u = hu, where e = g(u, u) — g(a -u,a-u) = a 2 - g(u, u) — a 2 ■ e, 
g(u) =g{a-u) =a-g(u). 

From the explicit form of h u 

- 1 

n = g ■ u®u 

e 

and the form of u under the transformation of r = r(r), u = a ■ u , a :^ 0, it 
follows that 

h u = g } - -a-u®a-u= (136) 

g(a ■ u, a ■ u) 

-%u = h 1L . 



g(u, u) 



Therefore, h u = h u . By that e = a 2 ■ e, e :— g(u,u). 

Remark. Since the vector field u = (1/a) ■ u is collinear to the vector field 
u it is obvious that the corresponding to u and u projective metrics should not 
change when u change to u. Nevertheless this conclusion should be proved. 



6 Perfect fluids 
6.1 Introduction 

1. A fluid for which its thermodynamic properties and its viscosity are con- 
sidered as not important and could be neglected is called a perfect fluid. This 
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means that the motion of a perfect fluid does not depend on heating and viscous 
processes The model of a perfect fluid appears as the simplest model of a 
fluid. A perfect fluid has no stress tensor (the stress tensor is equal to zero) and 
no orthogonal to its motion (its velocity) conductive energy flux density and 
conductive momentum density |12| . 

A perfect fluid is usually defined and characterized by two different ways: 

(a) By means of the Euler-Lagrange equations as equations for a fluid in 
static or stationary state (Euler's equation, Bernuli's equation). 

(b) By means of energy-momentum tensors with vanishing stress tensor, 
conductive momentum density, and conductive energy flux density. In the rel- 
ativistic continuous media mechanics the second way is preferred. The reason 
for that is that the Euler-Lagrange equations appear in general as sufficient but 
not as necessary conditions for the existence of an energy-momentum tensor of 
a given type. 

Let us now consider the characteristics of a perfect fluid more closely. 
Definition. A perfect fluid is a continuous media with an energy-momentum 
tensor G of the type 

G := (p G + i • k ■ L) ■ u®g{u) - L ■ Kr , (137) 

where 

Pg = - 2 -[g{u)](G){u) , 

k = -■[g{u)]{Kr){u) , 

e 

Kr — g l A ■ di ® dx^ , e = g(u, u) 

The scalar invariant L is a given Lagrangian invariant, interpreted in the 
case of fluid mechanics as the pressure p of a fluid, i.e. 

L:=p . (138) 

Proposition 2. The necessary and sufficient conditions for the existence of a 
perfect fluid are the conditions: 

G 7f = fe 7T = , 

G s = k s = , (139) 
G £ = k S = . 

Proof: 1. Necessity. From the general representation of an energy- momentum 
tensor G in the form [S], 0] 

G = (p G + - ■ k ■ L) ■ u® g{u) - L-Kr + 
e 

+u(E> g( k Tr) + fe s<g> g(u) + ( k S)g , 
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it follows that the condition 

u®g{ k Tr)+ k s ® g{u) + { k S)g =0 , (140) 

appears as a necessary condition for the existence of a perfect fluid (in accor- 
dance to the definition for a perfect fluid). 
After contraction of the last relation: 

(a) with g(u) from the left side, we obtain 

e • fc 7r = , k TT = for e ^ , 

because of g{u k s) = and g(u, { k S)g) = 0. 

(b) with u from the right side, it follows that 

e • k s = , fc s = for e ^ , 

because of g( k ir, u) = 0, ( k S)g(u) = 0. 

For k TT — and fc s = 0, it follows that k S = 0. Therefore, the necessary 
condition is equivalent to the conditions (|140|) . 

2. Sufficiency. If k n = 0, fe s = 0, and k S = then 

it (g) g( fe 7r) + fc s (g g{u) + ( k S)g = 
and the energy-momentum tensor G takes the form 

G = {p G + - - k- L) - u® g(u) — L ■ Kr . 
e 

The energy-momentum tensor G could also be written in the form 

G = k G - L ■ Kr . (141) 

The tensor k G is called viscous (viscosity) energy-momentum tensor and k S 
is called stress tensor. 

A comparison of the last expression with the expression for G of a perfect 
fluid leads to the form of k G as 

k G=(p G + --k-L)-u®g(u) . (142) 

After contraction from the left side with g(u) and from the right side with 
u, we obtain 

[g{u)]( k G)(u) = (p G + --k-L)-g(u,u)-g(u,u) = 

e 

= c 2 ■ (Pg ^ k- L) = e 2 ■ p kG , 

e 

[g(u)]( k G)(u)=p G -e 2 + k-e-L , 
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P* G = -^-lg(u)]( k G)(u)=p G + l--k-L , (143) 

£ = ! ■ (p* G - P G ) • ( 144 ) 

Therefore, every Lagrangian invariant L obeying the above condition could 
be used for description a perfect fluid. 

If we interpret the Lagrangian invariant L as the pressure p, i.e. if L = p, 
in a dynamic system the condition for L appears as a state equation for the 
pressure p of a perfect fluid 

p=^-.{p kG -p G ) = 1 ±-.[g(u)]( k G)(u)-yp G -e . (145) 

2. The energy- momentum tensor G for a perfect fluid could also be written 
in the form 

G = {p G + --k-L)-u® g(u) — L ■ Kr = 

e 

= p G ■ u<gi g(u) H — ■ L ■ {— ■ [g(u)](Kr)(u) ■ u <g> g(u) — Kr} (146) 
e e 

because of k = (1/e) • [g(u)](Kr)(u). 
Therefore, 

1 1 

G = p G -u®g(u) L-{Kr [g(u)](Kr)(u) ■ u <g> g(u)} , 

(G)g = p G ■ u <g> u - ^ • L ■ {{Kr)g - ^ • [g(u)](Kr)(u) ■ u <g> u} . (147) 
Special case: (L„, g)-space: S :— C, k = 1. 

(G)<? = (p G + ^ • L) • m® u- X • 5 (148) 

because of 

(A-r)ff - • g jfe • ft ® ft = 5 lfe • ft ® ft = 3 . (149) 
The tensor (G)g could be written in the form 

(G)g = {Pg + \- l ) ■ u®u- L -g = 

= p G ■ u® u — L ■ (g — - ■ u® u) = 

= p G -u®u- L-h u , (150) 

h u = g — — ■ w ® u . 

e 



Therefore, 



9[(G)g] = gjj ■ G ij = p G ■ g[u ®u]-L- 9 [h u ] 
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Since 

g[u <g> u] = g(u, u) = e , g[h u ]=n—l , (151) 
we have for g[(G)g] 

g[{G)g]=p G -e-{n-l)-L (152) 

and 

L=-^ J .(p G .e-g[(G)g}) . (153) 

If L = p, i.e. if L is interpreted as the pressure p then p could be expressed 
in the form 

P=^-(p G -e-g[(G)g}) . (154) 
For an energy-momentum tensor G with 

g[(G)g] = (155) 

we have 

Therefore, if a perfect fluid is described in a (L n , g)-space by means of the 
corresponding for this fluid energy-momentum tensor G, obeying the additional 
condition (having the additional property) 

g[(G)g] = 9ij ■ G" !ltj ■ G lj = 9ij ■ G ij = (157) 

then the pressure p of the perfect fluid is proportional to the rest mass density 
Pq of the corresponding energy-momentum tensor G. 

Special case: (L n , <y)-spaces, dimM :— n := 4, e := eo := const. :^ 0, 

g[{G)g] = o. 

P=T i -pG- e ■ ( 158 ) 

If the pressure p and the energy-momentum tensor G are given, we could 
find out the value l u = \ u \ = \ g(u,u) I 1 / 2 of the velocity vector of the perfect 
fluid using the form of e = ±1^ 

e = ±ll = — -{(n-l)-p + g[(G)g}} . (159) 
Pg 

Remark. The sign before b\ is determined by the signature Sgn g of the 
metric tensor g. If Sgn g > then e = — 1\. If Sgn g < then e = +/ 2 . 
Special case: (L n , g)-spaces, g[(G)g] :— 0. 

±ll = — -(n-l)-p , (160) 
Pg 



(n- !)■ I — I • (161) 
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For p > 0, p G > 0, we have 



l u = J(n-l)-— . (162) 
V Pg 

Special case: (L n , g)-spaces, g[(G)g] := 0, n := 4, p > 0, ,o G > 0. 

l u = V3- . [~^~— ■ (163) 
V Pg 

For a (L„, g)-space, since g[G(g)] = g[(G)g], we obtain for a perfect fluid, by 
the use of the explicit form of G, the relations: 

9[{G)g] = (pG + --k-L)-g[u®u]- P - g[(Kr)g] . 

e 

Since 

g[u®v\ — g(u,u) = e , 
g[(Kr)g] = g^ • g) • g* = g^ • <?* = • g mk • f m = 

n m . fJ _ fin _~f 
— yj J m — J m — J > 

it follows for g[(G)g] 

9[{G)g] = p G -e + k-p-p-J = 

= p G -e-(f-k)-p . (164) 

Then 

P==^—-(p G -e-g[(G)g}) , (165) 
f-k 

e= — .{(f-k)-p + g[(G)g]} . (166) 
Pg 

For (L n , g)-spaces k = 1 and / = n (dimM = n) . 

7 Lagrangian formalism for a perfect fluid in 
special types of (L n , g)-spaces and in (L n ,g)- 
spaces 

Let a Lagrangian density of the type 

L =V /I 5'Pfej.ftP; l E P !t .«') (167) 

with 

P -=Po+Pi(p) + ao ■ p-e+p 2 (up) (168) 
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be given, where 

P 2 ( u /°) = P2( uk P,k) > e = g(u,u) , (169) 

p = const. , a = const. , 

p = P {x k ) e <E>° (M) . 

On the basis of the method of Lagrangians with covariant derivatives (MLCD) 
the whole scheme of the Lagrangian theory corresponding to the pressure p as 
Lagrangian invariant could be found in its explicit form. 

7.1 Euler-Lagrange's equations 

7.1.1 Euler-Lagrange's equation for the invariant function p 

For the invariant function p the Euler-Lagrange equation is 

dp , dp . dp 
dp dp ti dp A 

where 

q = T ki k ---g kl ■ 9ki ;i + 9 l k ■ 9i-,i ■ 

By the use of the relations 

dp dp 1 
op op 

dp dp 2 (u k -p k ) _ dp 2 d k 

dp., " dp, d{up)' dp ;i {U ' Pk) 



dp 2 k A _ dp: 



d(up) ' U ' 9k 8(up) " ' 

^El + ao . e _ ( d ^ 2 . u l ) - + -q- ■ u l = 

dp d{up) ,% d{up) 1 

dpi . u dp 2 dp 2 t 



2 



dp 2 



■ u l , = 



d(up) 

the Euler-Lagrange equation for p follows in the form 
dP\ . r / dp 2 . dp 2 t dp 2 

Special case: p 2 := bo ■ up = bo ■ u % ■ p t i, bo — const. ^ 0. 

d P~2 _ K 
d{up) ° ' 

^ + ao ■ e + b ■ q t ■ u l - b ■ u l = . (172) 
dp 
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7.1.2 Euler-Lagrange's equations for the velocity vector u 

The Euler-Lagrange equations (ELEs) for the vector field u could be found in 
the form 

— = a ■—+ ^ 2 ■ ^ UP) = (173) 
du l du l d(up) du l 



Since 



de k 
d(up) 

~faT = p > i ' 

dp k , &P2 n 

the ELEs for u follow in the form 

2.a .p.g rk .u k + -^yp, i = . (174) 

Special case: p 2 := b ■ up = b ■ u % ■ p t i, 6 = const. ^ 0. 

2 • «o ■ p ■ 9m ■ uk + b o ■ P,i = , 

P, = -2-^-p-g- k -u k . (175) 

7.1.3 Euler-Lagrange's equations for the metric tensor g 

The ELEs for g could be written in the form 

d P , 1 _ jw 



dgu 2 



-•p-. 9 fe ' = 0. (176) 



From the relations 



dp de d 

Iki 



dgu dgu dg k - 



I 1 d 9ij 

= a • p-u-u 3 ■ , 
ogu 

dp 



dg, 



ki 



a ■ p ■ u k ■ u l 



we obtain 



dp , 1 „ M „ „,k i , 1 „ „fci 



+ • P • 5 = a ■ p ■ u K ■ u l + - ■ p ■ g Kl = 



%d 2 2 

S w = -2 • a • - • u fe • u ( , p :^ . (177) 
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7.1.4 Corollaries from the Euler-Lagrange equations for the metric 
tensor g 



From 



a • p ■ u k ■ u l + i • p ■ g kl = 



after contraction with g-^ we obtain 

1 2 
a ■ p ■ e + — ■ n ■ p = : p= • a ■ p ■ e . (178) 

On the other side, 

p = Po+Pi(p) + aa- P- e+p 2 (up) . 
From both the expression for p, the condition 
2 

(1 + -) • a ■ p ■ e +Pi(p) +P 2 ( U P) = ~Po = const., n>0 , 

follows. The metric tensor components g kl could be represented by means of 
the components h kl of the projective metric h u as 

g kl = h kl + --u k -u l 

e 

and on the other side as 

g kl = -2 • a • - • u k ■ u l . 
P 

From the last two relations, it follows the relation 

h kl = -(- + 2-a --)-u k -u l 
e p 

which contradicts to the property of h u to be orthogonal to the vector field u, 
i.e. the explicit form of h kl contradicts to the properties of h u 

h u [g{u)] = Q , h^-gj- k -u k = 0, 

because of 

h kl -9j m -u m = -(± + 2-a ^)- 9lm -u™-u k -u l = 

= -(i + 2-a --)-e- U fc ^0 , (179) 

e p 

e ? 0, - + 2-a -^^0 , 
e p 

Remark. If 

1 P 
- + 2 ■ o ■ - = 

e p 
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then 

-=~o . p=-2-a -e-p , h kl =0. 

p 2 ■ ao ■ e 

The expression (I179f) means that a projective metric h u , orthogonal to u, 
cannot exist if the ELEs for g are fulfilled. Therefore, we cannot consider the 
sub space, orthogonal to the vector field u, if the ELEs for g are valid. From 
this point of view, the metric tensor field g could not be considered as dynamic 
field variable. It should be assumed to be given as a non-dynamic field variable. 

7.2 Energy-momentum tensors 

7.2.1 Generalized canonical energy-momentum tensor 

?i S = Jj-/P,i-P-9i • (180) 



Since 



we have 



dp = dp 2 
dp d d{up) 

- dP2 n .„*_„.„./ 



•/<,•"' /-•//; ■ (181) 



d{up) 

Special case: p 2 := bo ■ up = bo ■ u 1 ■ p^, bo — const. ^ 0. 

6i j = b -p A -u j -p -gl . (182) 
If the ELEs for the vector field u are fulfilled then 

0i j =-2-ao-p-g m -u k -u j -p-g{ . (183) 

7.2.2 Symmetric energy-momentum tensor of Belinfante 

The symmetric energy-momentum tensor of Belinfante has the simple form 

-/.' i'-n! ■ (184) 

7.2.3 Variational energy-momentum tensor of Euler-Lagrange 

The variational energy-momentum tensor of Euler-Lagrange has the form 



i 3 = vQi 3 + gQi 3 , (185) 



7; j j 1 Sp k Op 
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ogyk og 3 _k 

= -2 • a ■ p ■ gjj: ■ u k ■ v? . 

By the use of the above relations, the variational energy-momentum tensor 
could be written as 

= 2 • a ■ p ■ ■ u k ■ u 1 + ■ P>i ■ u> - 

-2.a Q -p- grk -u k -u> = ^yP,i-ui . (187) 

Special case: p 2 := b$ ■ up = b$ ■ u 1 ■ p^, bo — const. ^ 0. 

Q,' b.-p,-,r . (188) 
If the ELEs for the vector field u are fulfilled then 

Q, ' 2-,/ , • /) . / y-.,/ / •//' . (189) 

7.3 Invariant projections of the energy-momentum ten- 
sors 

7.3.1 Invariant projections of the generalized energy-momentum ten- 
sor 6 

(a) Rest mass density p B . 

Pe = \ ■ Sk '' ■ uj ■ u k = 

1 dp~2 k 1 i i 1 k 

= - ■ — — - ■ u ■ p u ■ k ■ p , k — - ■ U-r ■ u , 

e d{up) e e k 

Pe = i.[^.^.p ifc -*.p] . (190) 

Special case: p 2 := b a ■ up = b a ■ u l ■ b — const. ^ 0. 

Pe = - e -[b -u k ■ Pk -k-p] . (191) 

If the ELEs for the vector field u are fulfilled then 

Pe = ~ ■ (2- a ■ p-gjrf u k ■ u l + k ■ p) . (192) 

Special case: (L n , g)-spaces: S = C : k = 1. If the ELEs for u are fulfilled 
then 

p e = -(2-a -p+-) . (193) 

e 
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(b) Conductive momentum density 6 tt. 

1 , 



- • k 0i k ■ u r • ft" = 

e K 



- 0(„„) ■ (194) 

Special case: p 2 := b ■ up = b ■ u l ■ 6 = const. 7^ 0. 

V = fe • PjJ • ^ . (195) 
Special case: (L n , g)-spaces: S = C : k = 1. 

n= WF)' Py ■ (196) 

If the ELEs for the vector field zt are fulfilled then 

V = -2-a -p- fe - U m -/i l4 = (197) 
= — 2 • a ■ p ■ uj ■ -h . 

Special case: (L n , g)-spaces: S = C : k = 1. If the ELEs for w are fulfilled 
then 

V = -2-ao- p-gi m -u m -h li = . (198) 

Therefore, in a (L„ , <?)-space the conductive momentum density 8 7r vanishes 
if the ELEs for the vector field u are fulfilled. 

(c) Conductive energy flux density 8 s. 

e r = l-h ij -gjj:- k 0i k -u l = (199) 

= • 

The conductive energy flux density e s is equal to zero for the given La- 
grangain density. 

(d) Stress tensor e S. 



e s^ = ti k - m - k e m l -h m i = 

- ®P 2 . o . J . h ik ■ a- ■ - 

~ d{u P ) p - m m 

' kh p !m .h^.h ik . m .u l = . (200) 



d (up) 

The stress tensor e S is equal to zero for the given Lagrangian density. 
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7.3.2 Invariant projections of the symmetric energy-momentum ten- 
sor of Belinfante S T 

(a) Rest mass density p T . 

= ---k-p . (201) 

e 

Special case: (L n , g)-spaces: S = C : k = 1. 

p T = ---p . (202) 
e 

(b) Conductive momentum density t tt. 

T ^ = \- sk Ti k -u^-h li = Q , sk Tt k = . (203) 

(c) Conductive energy flux density T s. 

T s l = - ■ fc« ■ m ■ sfc T ; k ■ u l = , sfc T z k = . (204) 
e J 

(d) Stress tensor T 5. 

T S^ = h ik ■ g Yl ■ sk T m l -h™i=0 , sk Tt k = . (205) 

7.3.3 Invariant projections of the variational energy-momentum ten- 
sor of Euler- Lagrange 

(a) Rest mass density pq. 

I>Q = -^•9iJ'" i '5t'-«* = -^-^y(>,<,-»'-«i-«' = 

= 4Hy ■ (206) 

Special case: p 2 := b ■ up = b n ■ u 1 ■ p^, b = const. ^ 0. 

p Q = -- e -b -p >k -u k . (207) 

If the ELEs for the vector field u arc fulfilled then 

1 



Pq = ~~ e ■ (-2 • a -p-gji- u l ) ■ u k = 
2 



a ■ p ■ -u l ■ u k . (208) 
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Special case: (L n , g)-spaces: S = C : k = 1. If the ELEs for u are fulfilled 
then 

p Q = 2 • a Q ■ p . (209) 

For a := \ the rest mass density pq is identical to the introduced invariant 
function p in the Lagrangian invariant L = p. 
(b) Conductive momentum density ^n. 

<V = -I. Q^.u h U = ^.^ prU k. u h U = 
e ' k e d(up) k 

- ^ -Pyh U ■ (210) 



d{up) 

Special case: p 2 := b ■ up = b n ■ u l ■ 6 = const. 7^ 0. 

Q ir l = -b - P r h H . (211) 

If the ELEs for the vector field w are fulfilled then 

Qjr i = 2-a -p-g lm -u m -h li • (212) 

Special case: (L n , g)-spaccs: S = C : k — 1. If the ELEs for w are fulfilled 
then 

<V = 2 • a • p • gi m ■ u m ■ h H = . (213) 
(c) Conductive energy flux density Q s. 

Q s l = -- ■ g-jk ■ Qi k ■ u J = 

e J 



p,-u k -W -9^-^ = . (214) 



(d) Stress tensor Q S. 



h ik -a ... • „' ■ 

. „ . hrnj . Life 



9(up) 



• p,™ • • • m ■ « ' = • (215) 



We can summarize the results for the projective quantities if the ELEs for 
the vector field u are fulfilled in a (L n , <?)-space and in a (L n , <?)-space in the 
following table 



Quantity 




e 7f 




e s 


Pt 


T— 


T— 

J s 


T s 


PQ 




Q s 


Space 






















(L„,5)-space 




^0 








^0 













^0 





(L n ,#)-space 


^0 











^0 











^0 
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It is obvious from the table that a perfect fluid described by the use of the 
given Lagrangian density L could exist only in a (L„,g)-space and, in general, 
it could not exist in a (L n , g)-space. In special cases of (L n , g)-spaces with 
conformal contraction operator S with S(dx l ,dj) — f l j := ip{x k ) ■ g* there are 
analogous relations as in (L n , g)-spaces. In a (L n , <?)-space: 

= {p e + --p)-u®g{u)-p-Kr , Pg = -(2 • a • p + ?■) ,(216) 

e e 

S T = -p-Kr , p T = -- , (217) 

e 

<3 = -PQ-u®g{u) , p Q = 2-a -p , (218) 

6* = -2- au ■ p- u® g(u) - p ■ Kr , (219) 
S T = -p-Kr , (220) 
Q = -2 • a • p- u®g(u) . (221) 

If we chose the arbitrary constant ao ^ as a := — \ then 

p e = p-- e , (222) 

Pt = -- , (223) 

e 

Pq = ~P ■ (224) 

The scalar invariant function p is proportional to the rest mass density pq 
of the variational energy-momentum tensor of Euler-Lagrange. The rest mass 
density p T of the symmetric energy-momentum tensor of Bclinfante is equal up 
to a sign to the rest mass density generated by the pressure p of the system. 
The rest mass density p e corresponding to the generalized canonical energy- 
momentum tensor is equal to the difference of the rest mass p and - (if ao = — \). 

The generalized canonical energy-momentum tensor could be considered as 
the sum of the symmetric energy-momentum tensor of Bclinfante and the vari- 
ational energy-momentum tensor of Euler-Lagrange 

8= S T + Q . (225) 

The last relation represents the second covariant Noether identity for a per- 
fect fluid. 

In the further considerations we will assume the existence of a perfect fluid 
in a (L n , g)-space or in its special cases. 



8 Navier-Stokes' and Euler's equations in (L n , g)- 
spaces 

8.1 Navier-Stokes' equations 

In a (L n , g)-spacc (S = C), if the Euler-Lagrange equations for the scalar func- 
tion p and for the vector field u are valid, we have the relations for a perfect 
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fluid in a co-ordinate basis: 



Since 



e 

+[(Pe + ~ ■ P),j ■ u3 + (Pe + ~ ■ P) ■ u3 ;j] ■ ul ~ 
-P,j-9 t3 + {Pe + --P)-9 a -gixk-v? -u k . (226) 



p = -(2-a o -p+^-p) , 



we have further 

g ik -0k j ;j = -2-a -p-a l - 

-2 ■ a • (pj • u J + p- v? .j) ■ u 1 - pj ■ g l] - 
-2-a -p-g il -g l0 , k -u j -u k . (227) 

The Navier-Stokes equation could be written in the forms 

(a) h u [g{56)] = , (228) 
(b) /'/,•//'"• = (229) 

Since (see the above consideration of Navier-Stokes' equation) 

p = p g + - -p= -2-a - p , fc = l , 5Kr = , (230) 

e 

the Navier-Stokes equation for a perfect fluid in a (L n , <?)-space could also be 
written in the forms [see the special case of (L n , g)-spaces with G W := k n = , 
G s:= k s = 0, G S := k S = 0.} 

-2 • a ■ p - a - g(Krp) + ^ • [2 • a ■ p • g(u, a) + g(u, g(Krp))] ■ u- 

-2-oo-p ■ h u [(V u g)(u)} = . (231) 

2 • a ■ p ■ a = -g(Krp) + i • [2 • a ■ p ■ g(u, a) + g(u, g(Krp))] ■ u - 

-2-oo-p -h u [(W u g)(u)] , (232) 

2 • a • p -a = -g(Krp) - 2 • a ■ p ■ h u [(V u g)(u)} + f ■ u , (233) 

where 

/ = - • [2 • a ■ p ■ g(u, a) + g{u, g(Krp))} . (234) 

e 

Special case: U n - and V^-spaces: V^.g := for V£ € T(M). 

2 ■ a ■ p -a = -g(Krp) + f ■ u . (235) 
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a = - 7r — • =2 • g(Krp) - - 2 ■ h u [(V u g) («)] . (238) 



If we apply now the general consideration for the transformation of the 
proper time r of the vector field u = d/dr we can find the Navier-Stokes equation 
in a new form by the use of the relations 

a — a 2 -a , a = a® ■ exp[f • dr] , ctQ = const., 

^ = 2 • a ■ p ' / ' 

a = ~ - 1 ■g(Krp)-h u [(V u g)(u)] + f-u , (236) 
Z ■ a ■ p 

a — f ■ u = a 2 - a , 

a 2 -a=--^—-g(Krp)-h u [(\7 u g)(u)} . (237) 
Z ■ a ■ p 

For the new proper time parameter r the Navier-Stokes equation has the 
form 

2-h- P -k 2 -~ 9{Krp) -h 

Since the metric tensor field ~g is an arbitrary given contravariant metric, 
corresponding to the covariant metric tensor g, we can introduce a metric g := 
a 2 ■ g, conformal to the metric g. Then the corresponding contravariant metric 
~g will have the form 

V=^2-9 • (239) 

a 

For the new introduced metric g the Navier-Stokes equation could be written 
in the form 

a = - - J— • f (Krp) -h u [^- (V u g) («)] . (240) 
Z ■ aa ■ p a 

On the other side, 

u = a-u , h u = h w , (241) 
V u g = V 3 . s .g = a • V^.g , (242) 
(V u g){u) = (Va-ug){a-u) = a 2 -(X7ug){u) , (243) 

i-(V u5 )(«) - (V^)(u) . (244) 
a 

By the use of the above relations the Navier-Stokes equation takes the form 

o = - ■ V(Krp) - h^KVug) (u)} , (245) 

Z ■ a ■ p 

which is exactly the generalization of the Euler equation for (L n , g)-spaces. In 
U n - and V^-spaces (V^g = 0) and for a := | • e with (e = =Fl) the equation 
takes its standard form 

a = e ■ - ■ Tj(Krp) (246) 
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or in a co-ordinate basis 

a' = s-yg ij -pj , e = Tl • (247) 

Remark. The sign of e is depending on the sign of the pressure p defined 
by different authors as extrovert (acting out) or introvert (acting into) quantity 
with respect to a dynamical system. 

9 Conclusion 

In this paper the Euler equation for a perfect fluid is considered as a special 
case of Navier-Stokes' equation in spaces with affinc connections and metrics. 
By the use of a special type of a Lagrangian density the Navier-Stokes equations 
are obtained. It turns out that these equations could be considered as equations 
for a perfect fluid in spaces with one affinc connection and a metric [(L n ,g)- 
spaces] and only in special cases of spaces with contravariant and covariant affinc 
connections and metrics [(L n , g)-spaces]. 

The Euler equation is obtained by the use of entirely unconstrained La- 
grangian formalism represented by the method of Lagrangians with covariant 
derivatives (MLCD). The additional conditions by this method are related to 
the transformation properties of the proper time of a mass element (particle) 
moving in space-time. They are not related to any constraints of the variational 
principle used for finding out the Euler-Lagrange equations for a perfect fluid. 
Therefore, we can consider unconstrained variational principles in hydrodynam- 
ics but under the condition for more careful use of the notion of proper time (or 
of time parameter) and its transformation properties. It seams that the notion 
of proper time (as the notion of time in general) could be much more important 
in physics than assumed until now. 
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